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We inquire about the possible coexistence of macroscopic and microstructured phases in random 
Q-block copolymers built of incompatible monomer types A and B with equal average concen- 
trations. In our microscopic model, one block comprises M identical monomers. The block-type 
sequence distribution is Markovian and characterized by the correlation A. Upon increasing the 
incompatibility \ (by decreasing temperature) in the disordered state, the known ordered phases 
form: for A > A c , two coexisting macroscopic A- and _B-rich phases, for A < A c , a microstruc- 
tured (lamellar) phase with wave number fc(A). In addition, we find a fourth region in the A-x 
plane where these three phases coexist, with different, non-Markovian sequence distributions (frac- 
tionation). Fractionation is revealed by our analytically derived multi-phase free energy, which 
explicitly accounts for the exchange of individual sequences between the coexisting phases. The 
three-phase region is reached, either, from the macroscopic phases, via a third lamellar phase that is 
rich in alternating sequences, or, starting from the lamellar state, via two additional homogeneous, 
homopolymer-enriched phases. These incipient phases emerge with zero volume fraction. The four 
regions of the phase diagram meet in a multicritical point (A c , Xa), at which A-B segregation van- 
ishes. The analytical method, which for the lamellar phase assumes weak segregation, thus proves 
reliable particularly in the vicinity of (A c ,Xc). For random triblock copolymers, Q = 3, we find the 
character of this point and the critical exponents to change substantially with the number M of 
monomers per block. The results for Q — 3 in the continuous-chain limit M — > oo are compared to 
numerical self-consistent field theory (SCFT), which is accurate at larger segregation. 

PACS numbers: 64.60.-i, 82.35.Jk, 64.60.De, 64.75.Va, 64.70.km, 64.60. Kw 



I. INTRODUCTION 

Random A-B block copolymer melts represent an in- 
teresting class of materials both for applications, due 
to their molecular self-organization for templating struc- 
tures on the nanoscale as well as for everyday materials, 
and theoretically, as mult i- component systems with com- 
peting interactions and a complex phase behavior [IH1]- 

For copolymer mixtures, phase separation was first ad- 
dressed by Scott [5] within a mean-field theory of multi- 
component demixing based on Flory-Huggins theory (see, 
e.g., [5]). Scott computed the limits of stability of the 
disordered, mixed state against macroscopic phase sep- 
aration for arbitrary distributions of chain composition 
(overall fraction of one monomer type). The coarse- 
grained description of ref. [5], which disregards the con- 
formations of individual chains, was subsequently ex- 
tended by Bauer [7] to assess the coexistence of multi- 
ple homogeneous phases and the equilibrium transition 
lines. The method was applied to random copolymers 
by Nesarikar et al. [8] , who computed the phase diagram 
for various chain lengths and average compositions. The 
system is treated as a multi-component mixture, with 
components distinguished solely by composition. Upon 
increasing the incompatibility, successive separations into 
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a growing number of homogeneous phases with different 
compositions are observed. 

Taking into account the internal structure of the 
chains and the block-type sequences in a melt of ran- 
dom block copolymers is crucial for the description of 
microstructured states (often termed microphase separa- 
tion) [2J [31 I5HTT] , see the example in the right panel of 
fig. [I] Fredrickson, Milner, and Leibler [5] formulated 
a microscopic model for random block copolymers, with 
one block composed of M identical monomers and with 
the block-type sequence distribution parameterized by a 
correlation A. Based on this model, they derived a mean- 
field free energy of Landau form in the limit of many 
blocks, Q —> oo. The resulting phase diagram shows an 
isotropic Lifshitz point, separating a line of instabilities 
with zero wave number (macroscopic phase separation) 
from a line of instabilities with finite wave number (mi- 
crophase separation) , cf. the lines in fig. [2] 

Several attempts have been made to go beyond mean- 
field theory and to consider the effects of fluctuations, 
predicted to be particularly important for the instability 
at finite wave number [TSJ, [T3j . Whereas the early works 
[THUS] deduced complete stability of the disordered state 
against microphase separation, it was later shown that 
proper inclusion of a local term in the Landau- Wilson free 
energy restored microphase separation [16] . The transi- 
tion was found to be weakly first order, yet wavelength 
and amplitude of the microstructured phases matched 
the mean-field predictions [3] rather well. 
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Monte Carlo simulations for symmetric random 
copolymer melts with different numbers Q of blocks per 
chain were performed by Houdayer and Miiller jTTJ, [TB] . 
In contrast to the mean- field calculations [3] , macroscopic 
phase separation was found only for small Q (in a A-range 
shrinking with increasing Q), and further increasing in- 
compatibility in the two coexisting homogeneous phases 
resulted in a remixed state. The latter was interpreted as 
the coexistence of three phases, two homogeneous ones 
and a third microstructured one with symmetric com- 
position, as predicted for random diblock copolymers 
(Q = 2), by simulation [TS] and self-consistent field the- 
ory (SCFT) [H]. For Q = 3, the simulations [H] pointed 
to a three-phase coexistence with fractionation accord- 
ing to sequences: While the two homogeneous phases 
displayed a higher content of homopolymers, copolymers 
accumulated in the microstructured phase. 

In this paper we aim at an analytical theory for three- 
phase coexistence due to sequence-specific fractionation: 
According to its internal structure, in particular the 
number of bonded A-B contacts, a sequence class, e.g. 
AAB I BBA, may have different concentrations in ho- 
mogeneous and structured phases. Our global copoly- 
mer distributions are symmetric in A/B content, which 
causes the A-rich and B-rich phases in a macroscopi- 
cally separated state to map onto each other by permu- 
tation of A and B. The distributions of these two phases, 
though different in composition, are not called fraction- 
ated, since they preserve the global concentration of a 
sequence class, e.g., of AAB / BBA. Their A excesses 
of opposite signs result from exchange of A- and £>-rich 
subspecies only within one sequence class. The A-rich 
phase, for instance, successively substitutes BBA chains 
with AAB chains, inversely the Z?-rich phase. In con- 
trast, we define sequence-specific fractionation to alter 
the sequence (class) concentrations in parts of the sys- 
tem such that microphase separation is favored in one 
part, while macrophase separation persists in the other. 

Our main results are the phase diagrams as a function 
of block correlation A and incompatibility x (see, e.g., 
fig. [3] below) showing a three-phase coexistence region of 
two homogeneous and one lamellar phase. Additional 
information concerns the volume fractions, the wave- 
lengths, and the sequence distributions of the fraction- 
ated states, as well as the behavior at the multicritical 
point. Some results provided by the analytical method 
have been briefly presented in ref . |21) . Coming from the 
macroscopically phase-separated state, a lamellar phase 
emerges with zero volume fraction (called shadow) and 
with finite amplitude; similarly, coming from the lamel- 
lar state, two additional homogeneous phases appear as 
shadows. The nature of the multicritical point, where 
four states of the system meet, depends on the number 
M of monomers per block: For M < 7, the wave number 
of the incipient lamellar phase vanishes continuously on 
approach to the multicritical point, and the segregation 
amplitude vanishes linearly. For M > 7 and particularly 
in the limit of continuous chains, the wave number re- 



mains finite, giving rise to metastable regions on both 
sides. In this case, the critical exponent for the segre- 
gation amplitude is 0.5. Detailed sequence-concentration 
diagrams of the coexisting phases show the partitioning 
according to their morphologies. Except for at the mul- 
ticritical point itself, the shadow phase emerges with a 
finite deviation from the global, A-defined distribution. 
A numerical SCFT study for continuous triblock copoly- 
mers covers larger segregation amplitudes, but yields a 
similar phase behavior. 

The paper is organized as follows: The microscopic 
model is introduced in sec. [XXJ. Free energies of macro- 
scopic and microstructured phase separation and the 
sequence-specific correlators are derived in sec. Ill In 
sec. IV we construct the free energy of a fractionated 



state and discuss the resulting phase diagrams in sec. [V] 
SCFT as a complementary approach is presented in 
sec. |VT Sequence fractionation is addressed in sec. VII 



A discussion of the methods is given in sec. |VIIl| followed 
by conclusions and an outlook in sec. |IX| 



II. MODEL 

A. Symmetric random block copolymers 

We consider an incompressible melt of j = 1 , . . . , N lin- 
ear, random A-B block copolymers in a volume V (fig. [l] 
shows triblocks). All chains have degree of polymeriza- 
tion L = QM, each of the Q blocks comprises M identical 
monomers, 



block 1 block 2 
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M 



M 



block Q 
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(1) 



Both types of monomers or segments are assumed to have 
the same statistical length, b. To formulate the effective 
repulsive interaction between monomers of different types 
(see section II B |, we introduce an A excess variable qj(s) 
for the type of segment s on polymer j, which takes the 
values +1 for A and —I for B. The type sequences of 
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FIG. 1. (Color online) Cartoon of triblock copolymers and 
lamellar phase separation. 

symmetric random block copolymers are generated by a 
Markovian polymerization process with average A excess 
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q = 2p — 1 = (p is the global concentration of type-A 
segments) and block-type correlation 



A:= (l-2p AB )e [-1.+1] 



(2) 



of adjacent blocks along a chain [3]. Here, pab denotes 
the probability that a block of type A is attached to one 
of type B in the synthesis. Assuming homogeneity, p and 
A are independent of the position on a chain. Positive A 
signal a preference for homopolymers, and A = de- 
scribes ideal (uncorrelated) block sequences. This model 
for the synthesis amounts to choosing the simplest non- 
trivial distribution with only two parameters, p and A. 
With p = 1/2 in our case, the corresponding transition 
matrix M for the probability vector [pa{P), 1 — p A {P)] T 
(probabilities to find A, respectively B, at block jS) reads 



M = 



l+A 1-A 

2 2 

1-A l+A 

2 2 



(3) 



Its diagonalized form is used to compute the probabil- 
ities of individual sequences in the A-distribution, and 
moments of the A excess distribution. Once generated, 
the block sequences remain fixed, i.e., thermal averaging 
affects only the chains' center of mass locations and their 
conformations. For a finite number of different sequences, 
a concentration for each sequence is well-defined in the 
thermodynamic limit. Hence, for finite Q, the quenched 
disorder due to the fixed block types on one chain can be 
effectively translated to a multi-component system. 

A chain can contain to Q blocks of type A, which 
defines Q + 1 classes of chains. This classification in 
the "crushed polymer approximation" (see, e.g., [7], [8]) 
is sufficient to study the separation into homogeneous 
phases (see sec. HID below). However, it neglects dif- 
ferences in the sequence of the blocks, i.e., the spatial 
structure of the chains (for example, the average A ex- 
cess for both AAB and ABA chains is 1/3). The spatial 
correlation of types along a chain is essential for the for- 
mation of structured phases with nonzero wave numbers. 



B. Potentials 



The Hamiltonian H consists of three parts, 



(4) 



which reflect intra- and interchain interactions of the 
monomers on a mesoscopic level. Explicitly, for the for- 
mer we consider the connectivity of Gaussian or ideal 
polymers [22] acting between monomers on the same 
chain, for the latter excluded volume and incompatibility 



acting between all monomers, in units of k^,T: 

N L-l 



3=1 s=l 

N L 



2p, 



V E E u(\r h ( Sl )-r j2 (s 2 ) 

r ■ 1 „ _ 1 



(5a) 
(5b) 
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31,2,81,1 



(5c) 



where the primed sum in eq. ( 5c I is shorthand for the con 



strained sum in eq. ( |5b[ ) (the constraint can be dropped 
in the thermodynamic limit). Spatial variables r are di- 
mensionless, rescaled from physical positions R via 



2dR a 



1,2, 



(6) 



with b the rms end-to-end distance of a Kuhn statistical 
segment and d the spatial dimension. Accordingly, the 
constant dimensionless monomer number density is 



Po ■-- 



NL _ NLb d 
~V '~ V{2d) d / 2 ' 



(7) 



One effective segment of our model usually represents 
many physical monomeric repeat units, as to fulfill the 
prerequisite of statistical independence of subsequent 
bond vectors in the coarse-grained Gaussian chain model. 

The excluded volume interaction eq. ( 5b ) must be ac- 
counted for, even if we later perform the incompress- 
ible limit, since A excess and total density fluctuations 
are coupled. The pair potentials U(r), W(r) are sup- 
posed to be short-ranged, and we approximate them 
by 6 functions, neglecting short-wavelength fluctuations. 
Conceptually, Gaussian chain connectivity and compress- 
ibility are effective potentials, which are obtained after 
intcgrating-out microscopic degrees of freedom. They are 
chiefly of entropic origin and thus originally proportional 
to k^T. The Flory parameter x expresses in an empiri- 
cal way the local free-energy change per monomer due to 
A-B contacts compared to a surrounding of monomers 
of the same type with larger attractive potentials [6] . Its 
main part is usually enthalpic, such that in the normal- 
ized eqs. ([5]), x is inversely proportional to temperature, 
X oc T~ 1 , and increasing incompatibility is equivalent to 
cooling. In the following, U-qT is set to unity. 



C. Order parameter 

A convenient order parameter that detects separation 
into A- and B-rich domains (phases) is the thermal av- 
erage of the local excess of A segments [3] , 

N L 

°{r) = EE*( S ) 5 ( r ~ r j( s )) = M r ) ~ Qb{t), (8) 

j=l s=l 
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i.e., the difference of segment densities due to A and B. 
As a second field, we introduce the total segment density 

N L 

Q( r ) = EE^ r - r ^ s )) = Sa{t) + g B {r). (9) 

3=1 »=1 

With these fields, and in the limit W(r) — > 5(r), the 
incompatibility ( 5c ) takes the standard form |23) 



H x - 



,l'V ( <*(r) - ( te (r)-«£> 



X 

4f?o 



\ d r (a(r)y 



(10) 



Note that as a zero of the incompatibility energy we have 
chosen the homogeneously mixed state where the local 
densities of A and B coincide with their global fractions 
throughout the system. Analogously, the excluded vol- 
ume interaction eq. (5b I in the limit U(r) — > 8(r) is 



2Qo 



i d r(g(r))' 



(11) 



III. FREE ENERGY 



In order to assess the phase diagrams, particularly 
phase coexistence for random block copolymer melts, we 
compute the free energies of basic phase-separated states. 
The two important control parameters are the incompat- 
ibility x and the block-type correlation A. Figure [2] shows 
the topology of the phase diagrams we will derive below. 
As discussed already by Leibler and co-workers [2 , 
the disordered state of the symmetric melt becomes un- 
stable toward either macroscopic or lamellar phase sep- 
aration, depending on A. Between these two well-known 
states, a new state will be shown to become stable, viz. 
the coexistence of three phases: an ^4-rich one, a B-rich 
one, both homogeneous, and a lamellar phase. Coming 
from the macroscopically phase-separated state, the new 
phase is created by expulsion of chains with many A- 
B contacts from the homogenous cloud phases (for the 
terms "cloud" and "shadow" phase, see [35]) into a sub- 
system, which then displays deviations from the global 
A-distribution. This fractionation increases the A excess 
amplitude in the homogeneous phases. More explicitly, 
lamellae can appear in the new phase because the altered 
sequence distribution with fewer homopolymers gives rise 
to a maximum of the structure factor at nonzero wave 
number, whereas the structure factor of the global distri- 
bution favors macroscopic phase separation. Conversely, 
starting from the lamellar phase, homopolymer chains 
arc expelled into two new homogeneous phases. Thereby, 
for values of A, at which the global sequence distribution 
favors lamellae, homogeneous phases become stable in a 
subsystem, resulting in a fractionated state. 

First, we present the free-energy densities of homoge- 
neous phases and of lamellae separately. The expressions 
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FIG. 2. (Color online) Qualitative phase diagram of random 
block copolymers. Global instabilities of disordered melt: 
solid (red) line: macroscopic phase separation for A > A c , 
dashed (blue) line: lamellar phase for A < A c . Fractionation 
creates an in-between state with three coexisting phases. 



are deliberately kept as simple as possible to focus on the 
effect of varying sequence concentrations. In the follow- 
ing section, we go on to set up a fractionated multi-phase 
free energy, allowing for sequence distributions different 
from the global one, and discuss three-phase coexistence. 



A. Free energy functional 

Starting from the Hamiltonian of eq. Q, we aim at 
computing the canonical partition function 



Z = J Y[d rj ( s )c 



-n 



(12) 



Pair interactions are formally decoupled via functional 
integrations over the collective density fields a and g, 
and over two conjugated interaction fields a and g (with 
Fourier modes {07-, gk}) that restrict a to the A excess 
and g to the total segment density (cf. eqs. ^ and ([9])): 



Z = J V [a, a, g, g 

u 



(13) 



X « 



4NL 



2NL 



iakT-k + igkg-k z u [a, g] 



In this expression, the inner conformational integrations 
have factorized into single-chain partition functions Zj. 
All N u chains with a given block- type sequence v, which 
is characterized by the segment types {q l ,(s)}, contribute 
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the single-sequence partition function 

z v [a, g] := (14) 

(exp{-i^^(a fe9l/ ( S ) + g k )c-* kr ^\\ . 

\ ^ fc^O s ' I 

Here, ((. . .)) denotes the conformational average 
/ Vr(s)(. . .) exp {-i Efji 1 (r(s + 1) - r(.s)) 2 } 



and the single-sequence partition functions 



JVr(s) exp {-iE s L =xH« + !) " r ( s )) 2 } 



(15) 



for one Gaussian chain (cf. eq. ( |5a[ )). Combinatorial pref- 
actors 1/N V \, homogeneous contributions (fe = 0), and 
the conformational partition functions of noninteracting 
Gaussian chains have been divided out in eq. (13), since 



we are interested in the free energy of a global ordered 
state relative to the disordered, homogeneous state. 

In order to perform the saddle-point approximation, 
we choose to first integrate out the amplitudes of the 
physical fields in favor of the conjugated ones, contrasting 
with the procedere in, e.g., refs. (see the note below 

eq. (p3|). From eq. (13), we obtain the linear relations 



2NL „ NL „ 
<7fe = -i <7fc and g k = I Qk, (16) 

X « 

and, for convenience, rescale the conjugated fields as 

(Th := i — — — fu and p k := ui k (17) 

fe 2NL h A^L y ' 

before insertion into Z. The resulting partition function 
in saddle-point approximation is 

Z « jV[f,(j] exp|-iV/i[f,w]|, (18) 

with the effective Hamiltonian (per chain) 
h[r,iu}= (19) 
( XTkt- h - l&k&-k ) - ^PuhiZv [f,w], 



fc^O 



Zu T, OJ = 



(20) 



-ik-r(s) 



The probabilities p„ : = N v /N define, in the thermody- 
namic limit, the sequence distribution over the up to 
possible realizations of a random, binary Q-block copoly- 
mer. (For Q > 2, the actual number of different se- 
quences is smaller due to the symmetry with respect to 
the two ends, see below for triblock copolymers.) 

Anticipating small field amplitudes, the next step is 
to expand the effective Hamiltonian h eq. ( 19 ) into a 



series in both fields. Restricting ourselves to systems with 
global A-B symmetry, the expansion contains no terms 
of odd order in f (the field conjugated to the A excess), 
since for n odd, moments of the A excess distribution, 



E 

Si,....S n 



qv(si)...q v (s n ),ne®, (21) 



are zero. A sufficiently large compression modulus n will 
prevent instabilities with respect to fluctuations of the 
total density. Hence, we can eliminate their conjugated 
amplitudes uj k perturbatively in favor of T k and obtain, 
to lowest order, a quadratic dependence 



\' 2 



(22) 



^(fei.fc-fei) 



8NL E , 1L + D(L,k 2 ) 
fci^{o,fe} k v i ; 



f kl T k - kl +0(\f k \ A ) 



(see appendix [A] for conformational averages of exponen- 
tials and appendix b] for the correlators and D). 
Substituting back this relation, and in the incompress- 
ible limit, k — y oo , the consistent expansion up to fourth 
order in f yields the free-energy functional per chain, 



X 4 W fSWffc^^Mftrfci-^-^) 



384(A^L) 4 

fci,fe 2 ,fe3 
fel+fe2+fe3^0 

A 



£ (3 g(L>(fcl + fca)2) S^( kl ,k 2 ,k 3 )jr kl r k2 r k3 r. hl . k2 . h3 



Y! {s^Hkl k 2 2 ) - S{kl)S(k 2 2 )) f fel f_ fcl f fc2 f_ fc2 + O (|f fe | 6 ) 
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with J2k(- ■ •) := Sfc^o- ^he global second-order corre- 
lator S(k 2 ), called structure factor in the following and 
discussed in sec. |IIIB[ is given by 



S{k 2 ) = Y J Pu £ quisMs^e-^'^ (24) 

V S1,S2=1 



written as an average over intra-chain correlators S l/ (k 2 ) 
of single block-type sequences. These and the correlators 
S<< a \ and SfW are defined in appendix |bJ In our 

global sequence distribution, the probabilities p v will be 
confined to A-defined values (see eqs. ( p9| below), but 
can take arbitrary values in a fractionated subsystem. 
As suggested by the functional eq. ( 23 ) , we assign the 



conjugated field f the role of the order parameter, since 
at the saddle point level, to which we adhere, averages of 
f and a are identical (cf. eqs. (16) and (17)). However, 



correlations of the conjugated field are not proportional 
to those of the field itself, cf., e.g., [2B]. Therefore the 



vertices in eq. ( 23 1 differ from those of the functional of a 



in refs. [31 [53] (apart from differences due to restrictions, 
e.g., to continuous chains with many blocks, which we 
do not impose). For instance, second moments of the 
amplitudes of a can be recovered from those of f via 



(f kl T-k 2 ) n - {Tk^u (r-h 2 ) n ~ 



2NL 

X 



(25) 



where is the canonical average, eqs. (13), respec- 
tively eq. (18). 



Aiming first at the simplest description, and in the 
spirit of a Landau free energy, we ignore the wave- vector 
dependence of the fourth-order coefficients in eq. ( 23 1 , 



i.e., we evaluate the correlators in the limit fc r — > 0, (in 
sees. 



fo [r] = 



III C 2 and V C we will relax this approximation) 
S(k 2 ) 



X 

AN 2 L 



£ 1 



2L 



(26) 



X 



128iV 4 



ki,k 2 ,k 3 



Tfei 7~fa 2 7~hz 7~— h%— k2— fea 



+ (m 4 - mf) £ r fcl f_ fel f fe2 f_ 



fei ,fc 2 



with the moments m2, TO4 from eq. (21 ) 



B. Structure factor and multicritical point 

The second-order structure factor S (fc 2 ) for a distri- 
bution of sequences sets the limits of stability of the 
homogeneously mixed melt. For our global Markovian 



distributions, solely the correlation parameter A decides 
whether the maximum position of S (fc 2 ) is located at 
zero or at finite wave number. In the former case, the 
disordered state becomes unstable with respect to macro- 
scopic phase separation, in the latter case to microphase 
separation [2 a . Upon decreasing A, the maximum posi- 
tion of S (fc 2 ) becomes nonzero at a critical correlation 
A C (M), depending on the number M of segments per 
block. The corresponding point in the X-\ plane where 
the lines of macroscopic, respectively lamellar phase sep- 
arations meet is termed a multicritical point, since also 
the transition lines to three-phase coexistence must end 
here. 

For a A-distribution of Q-blocks with finite M, the 
global S (fc 2 ) can be calculated from the probabilities of 
all type combinations of two segments with a given in- 
trachain di stan ce (in blocks) using the transition matrix 
M (cf. eq. ( B6 ) in appendix [Bj) : 



S^fc 2 ) =QD (M,fc 2 ) 

2Ae~ Mfc2 sinh 2 (^f 
+ (1-Ac- Mfc2 )sinh 2 (f) 



(27) 



1 - 



(Ae- Mfe2 )' 
1 - Ac- Mfc2 



with the dimensionless wave number fc 2 := 6 2 fc 2 / (2d) and 
fc the physical wave numb er. The discrete Debye function 
D (L, fc 2 ) is given in eq. (Bl). 

In the following, we restrict ourselves to the case of 
symmetric random triblock copolymers, Q = 3. This 
system features six different species, which we group into 
only three different (classes of) sequences, 



homopolymers: LLL; 
copolymers: KLL; 

LKL, 



(28a) 
(28b) 
(28c) 



K,L £ {A, B}, K ^ L, according to unfavorable intra- 
chain A-B contacts. Generally, pairs of species like AAB 
and BBA are related by blockwise A-B permutation and 
have the same topology of intrachain A-B contacts, and 
thus the same structure factor. To label these sequences, 
the index 1 is assigned to homopolymer chains ( 28a ) , 2 to 



copolymer chains with two adjacent blocks of the same 
type ( |28b[ ) , and 3 to strictly alternating chains ( 28c ) . For 
a A-distribution, the sequence (class) concentrations are 



Pi (A) = 
Pa (A) = 
Ps(A) = 



(1 + A) 2 
4 ' 
1-A 2 
2 ' 
(1-A) 2 



(29a) 
(29b) 
(29c) 



At a critical correlation X C (M), we find the following 
transition from macroscopic to lamellar phase separation: 

a) For M < 6, the maximum position of S (fc 2 ) is at 
fcp = for all A > A C (M)) and grows continuously 
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from ko = when A falls below X C (M) (see fig. 
low) . The critical value of the correlation, A, 
is reached when the second derivative of S (k 2 
k = changes sign: 



§be- 



at 



A C (M) = - 



1 - 



1 

M 



M < 6 



(30) 



b) For M > 6, however, a second maximum of S (fc 2 ) 
at k > evolves already for A > X C (M) (see 
fig. [7]). Now, the critical value A c is the one at which 
the second maximum (associated with a metastable 
lamellar phase) attains a higher value than the one 
at k = 0, and is accessible numerically only. 



For continuous Gaussian triblocks (segments indexed 
by a contour parameter instead of an integer) with un- 
altered coil diameter, the structure factors are computed 
in the combined limit M — >• 00, b 2 —> 0, Mb 2 = const, 
abbreviated as limM-s-ooi preserving the finite number of 
blocks, here Q = 3, and the rms end-to-end distance 
-Rbiock = \Mb. In this case, the wave number is con- 
veniently rescaled with i?biock- For a A-distribution of 
continuous triblocks, the global structure factor is 



s (fc 2 ) := lim S (k 2 /M) IM 2 = 

2Ae- fc2 sinh 2 (tf) 
3gD(1 ' fc2)+ (l-Ac-^/4 



(31) 



1 - f Ae~ fc2 



1 - Ac- fc2 



each triblock sequence (class) defined in eq. ( 28 ) 
S^k 2 ) = D(3M,k 2 ) 



(33a) 



3M(1 



) 2e- k (1 - e~ 6Mk ) 



(1 



S 2 (k 2 ) = -D{3M,k 2 ) (33b) 
+ 2 (D(2M, k 2 ) + D(M, k 2 )) , 

^(fc 2 ) = D(3M,k 2 ) 

- 4D(2M,k 2 ) + 8D(M,k 2 ). (33c) 

While the maximum of Si(fc 2 ) is located at k = 0, the 
maximum positions of ^(fc 2 ) and S3 (k 2 ) at k > are 
due to the finite type-position correlation length within 
a chain of the respective sequence. 

The continuous-chain version of the homopolymer 
structure factor eq. ( |33a[ ) is 

Sl (fc 2 ):= lim S, (k 2 /M)/M 2 = ffD (3,fc 2 ), (34) 

again with k 2 := R 2 lock k 2 /(2d); similar expressions hold 
for s 2 (k 2 ) and S3(k 2 ). In the following, S v {k 2 ) or S (fc 2 ) 
refer to the discrete structure factors, and the number of 
segments, M, is usually not listed as an argument sepa- 
rately. The continuous versions are denoted with s^(k 2 ), 
s (k 2 ) etc. Since the number of sequences grows exponen- 
tially with Q, the explicit calculation of sequence-specific 
structure factors is practically limited to a comparatively 
small number of different sequence classes, i.e., to a small 
number Q of blocks per chain. 



C. Lamellar phase separation 



now with k 2 := i? 2 lock fc 2 /(2(i), and with the continuous 
Debye function 



9o(i,k 2 



l + £k 2 



fc 4 /2 



(32) 



Continuous triblocks realize case [fe)J consistent with the 
case of triblocks with M > 6 discrete segments. The 
wave number of the global ordered (lamellar) state, fc (A), 
jumps dis ■continuously to zero as A approaches A c = 
—0.464 from below. The lamellar phase persists as a 
metastable state for A > A c , as well as macroscopic phase 
separation for A < A c . Remarkably, we discover this dis- 
continuity of the global wave number for the broader class 
of triblock copolymers with M > 6 segments per block, 
whereas the literature on copolymer mixtures seems to 
report only the behavior [a)] (see, e.g., [201 HZ]), associ- 
ated with a Lifshitz point [28] . 

Since we need to address sequence distributions differ- 
ent from the A-defined one in the next section, we calcu- 



late the second-order structure factors from eq. (24) for 



In order to derive the free energy due to microphase 
separation, we insert for our order-parameter field f the 
simplest single-harmonic ansatz [24] : lamellae with wave 
vector fc i : — l^ol > and an amplitude ffc 



Tk = Tfc (S, 



(35) 



More than one single wave vector is not considered here, 
since the instabilities of the disordered state of symmet- 
ric copolymers are known to be toward homogeneous or 
lamellar phases. In the latter case, we additionally as- 
sume that A-B separation is weak. 



1. Simplified lamellar free energy 

Insertion of the above ansatz into the simplified func- 
tional eq. (26) yields the free energy of a lamellar phase, 



fo(ko,T ko ) 

= ^-( 

2iV 2 L 



(36) 



2L 



647V 4 
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which is valid only for incompatibilities L\ exceeding 

2L 2 



(37) 



the onset incompatibility. (As usual, we shall use L\ 
instead of x as one parameter of the phase diagrams, 
due to the scaling of the onset incompatibilities with L.) 

Minimization of the function eq. ( 36 1 with respect to 
the order-parameter amplitude gives 



16V 2 




Lx 2 (m 



(38) 



Variation with respect to the wave number of the insta- 
bility shows that the optimal ko is the maximum position 
of S(k 2 ). With the single-harmonic approximation of the 
profile, the lamellar free energy at Lx > Ixm(fco) is 



fn 



\ L 2 



2 

Lx 



7714 



argmax S(k 

k>0 



(39) 



The first two phase diagrams in sec. [V] are based on this 
simplified version of the lamellar free energy. 



2. Lamellar free energy with restored wave-number 
dependence of fourth order coefficients 



With this ansatz, the free energy of Landau form becomes 



f - LX 
/h <° " H 



^rJ Th+ ^92^r Th ' (42) 



which provides a good description of macroscopic phase 
separation for small values close to the continuous 
transition from the disordered state. 

However, the transition we aim at, from the macro- 
scopically phase-separated to a three-phase state, may 
occur at a value Lx considerably larger than the onset in- 
compatibility of macroscopic phase separation; see fig. [3] 
Thus instead of the free energy eq. ( 42 1 that relies on an 
expansion in fh, we prefer and are able to derive a closed 
expression (cf. appendix [C]) by ignoring the copolymers' 
internal structure, consistent with uniform mean fields. 
For random triblock copolymers, the free energy is 

fh =—rir ~ i 1 ~ Pi) In cosh — p x In cosh — , 



provided Lx > Lxh 



2 

7772 



18 



l + 8 Pl 



(43) 



with the homopolymer concentration, p\ — (1 — p2 — P3) 
(the indices 2 and 3 refer to the sequence classification 
eq. |28| ) needed in the description of a lamellar phase). 
Here, the amplitude fh is determined by 



I-Pl 

3 



tanh 



LxTh 
Qg 



Pi tanh 



Lxh. 

2go ' 



(44) 



Restoring the fc-dependence of the fourth-order terms 



of eq. ( 23 1 , and optimizing the amplitude at a given wave 
number k m , we arrive at the free-energy function 



fm,k n 



(40) 



-L 



4 ( S(km) 2 



V L 2 



L X 



(5 (a) (fc m ,fc m )) c(P) (h u u \<q(f)(k 2 k 2 ) 
D (L 4fc 2 ) rn'" 1 "' K m ) -\- a \ y n m ,n m ) 

given x > 2L/S(k^), Now, minimization with respect 
to fc m results in a wave number that additionally depends 
on the incompatibility, fco = &o(x)- 



2. Homogeneous multi-phase coexistence 

Within multi-component theory, the two homoge- 
neous, A- and _B-rich phases in a random triblock melt 
(formed at Lxh — 6 for A = 0) are followed by four 
homogeneous phases at higher incompatibilities (e.g., 
Lxh, 4 ~ 16 for A = 0). More than four phases are im- 
possible within this theory, since for Q — 3 there are only 
four different chain compositions (A contents). The fact 
that no triblock sequence is symmetric in A/B content 
might explain why, starting from the A- and the B-rich 
phase, a third, homogeneous phase balanced in A/B does 
not become stable. 



D. Macroscopic phase separation 



1. Coexistence of two homogeneous phases 



Macroscopic phase separation can be assessed with a 



real-space version of the free-energy functional eq. ( 26 1 . 
Accounting for the symmetry, the appropriate ansatz is 
for two phases with uniform fields f of opposite signs in 
equally sized regions Vh,i and Vh.2 of the system: 



f(x) 



-f h , x e Vh.,2 



,1^,11 = 1^,21 



V 

2" 



(41) 



IV. FRACTIONATED THREE-PHASE 
COEXISTENCE 

In the following, we show that both a global macrosopic 
and a global lamellar phase separation become unsta- 
ble toward three-phase coexistence due to fractiona- 
tion. In the former case, mainly alternating sequences 
are expelled from the macroscopically phase-separated 
state (cloud) to allow for a third lamellar shadow phase, 
whereas in the latter case mainly homopolymers are ex- 
pelled from the lamellar state (cloud) to allow for two 
additional homogeneous shadow phases. 
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A. Fractionation from two macroscopic phases 

Here, we start at block-type correlations A > A c and 
incompatibilities L\ > Lxh, i-e., from a global, macro- 
scopically phase-separated state comprising two homo- 
geneous A-rich, respectively B-rich, phases. At further 
increase of L\, a third, lamellar phase with zero average 
A excess will be created by fractionation: Predominantly 
alternating sequences (y — 2,3) with few homopolymers 
will remix in a volume fraction := /V of the 
system. For our symmetric distributions, the two homo- 
geneous phases coincide in the volume fractions, in the 
field values up to the sign, in the sequence (class) concen- 
trations defined in eq. ( 28 ) , and thus in the free-energy 



densities. Hence we can treat them as one effective state, 
and study their joint sequence exchange with a lamellar 
phase. 

The first term of the free energy relative to the homo- 
geneous, two-phase state, is written as a weighted sum of 
the free-energy densities /m of the conjectured lamellar 
phase with volume fraction v^ 2 \ and f^ 1 of the two ho- 
mogeneous, homopolymer-rich phases with joint volume 
fraction v^ 1 ' = 1 — v^: 

= ({«?>}) +(l-t,W)/f> ({n«} 



fs\ 

'(45) 

Here, the sequence concentrations in state (phase) P are 
( p) 

denoted as n\, , v — 1,2,3, P = 1,2. The free-energy 
densities of global ordered states alone (for which com- 
binatorial terms due to the sequence distribution can- 
cel) cannot completely describe the coexistence of differ- 
ent states that interact via sequence exchange. Hence 
there are additional entropic coupling terms: First, con- 
finement of the chains to the volume fractions of phase- 
separated subsystems gives rise to a loss 



A/, 



vol. red. 



_ w (2) mu (2) _ _ w (2)) ln (! _ „(2)) ( 46 ) 



of translational entropy compared to the global state. 
Second, the sequence-selective exchange between the two 
phases effects a combinatorial gain A/ com b. due to the 
possibilities to choose chains of each sequence in one sub- 
system out of the total, A-defined number Np u (X) (the 
factorials are approximated by Stirling's formula): 



A/ com b. 



j^n^ln 




_ P„(A) _ 





(47) 



(1 



- j n 



(1) 



With the above contributions, the free energy of the frac- 
tionated phase coexistence is 



/frac — /sum ~t~ A/ vo i red. ~l~ A/ com b. 



(48) 



Incompressibility and the global A-defined sequence 
distribution reduce the number of variables, given by the 



volume fraction z/ 2 ' of the lamellar phase and the con- 
( p) 

centrations n v , v = 1,2,3, P = 1,2: The homopolymer 

(k) 

concentrations n\ can be eliminated by the constraints 



i/=i 



,(p) = 



= 1 for each phase P = 1,2. 



Likewise, the concentrations riv in the homogeneous 
phases can be explicitly expressed in terms of the volume 
fraction and the concentrations in the lamellar phase via 
the constraint of global A-defined concentrations, 

t/ 2) 4 2) + (1 - v (2) )4 1} = PuW, v = 2, 3. (50) 
Thus left with three independent variables, we choose 



(49) 



them as v^ 2 \ 



,(2) 



and 



/2) 



for the purpose of studying 



fractionation starting from two homogeneous phases. 

Obviously, at a given block-type correlation A and a 
give n incompatibility Lx, the fractionation ansatz eq. 
(48) is reasonable only for values of the variables v^ 2 \ 

n2 > n< 3^ f° r which /f rac reaches lower values than the 
free-energy density /h of the global state: 



A/6, 



(2) (2) (2)\ . 



/frac — fh < 0. 



(51) 



For each set (A, x), the free-energy change A/ frac has 



to be minimized with respect to v 



(2) r,.?) r,.i 2 ) 



the region limited by eq. (51). To avoid overloading the 



ng ' within 



presentation, the functional dependence on A and L\ is 
suppressed in /f rac , as well as in the free-energy densities 
of the global homogeneous and lamellar phases. 

To obtain the free energy of the lamellar state with 
fractionation, we compute the structure factor eq. (24) 



and the moments eq. (21) with modified sequence con- 
centrations P2 — > n\ and p 3 —> n y 3 ' (pi = 1 — p 2 — ps), 

(2) 

which are then the explicit arguments of / rrl . Similarly, 
the free-energy density of the macroscopically phase- 
separated state with fractionation is computed w ith mod- 



(!) 



ified concentrations n L 
becomes a function of v^ 2 \ 



, such that, via eq. (50), f h 
''\ and 



(i) 



B. Fractionation from a global lamellar phase 

The assumed boundary curve between the three-phase 
coexistence region and one lamellar phase comprising the 
total system, is in our approach restricted to the re- 
gion A < A c of the A-x space. To access this region, 
a fractionation ansatz has to start from lamellae in the 
A-distribution, which tend to expel homopolymers on in- 
creasing x, a mechanism which will give rise to homo- 
geneous A- and i?-rich shadow phases. The fractiona- 



tion free energy is formulated in analogy to eq. ( 48 1 in 
terms of the free-energy densities of one effective homo- 
geneous shadow phase, in a volume fraction z/ 1 ), and a 
lamellar cloud phase, in a volume fraction 1 — i/ 1 ). For 



10 



> 0, both states attain sequence concentrations de- 
viating from the A-defined ones. The first part of the 



fractionation free energy corresponding to eq. ( 45 ) is 



fs 



(52) 



(i) 



Again, the constraints eqs. (49) and (50 1 of incompress- 



ibility and fixed global sequence distribution reduce the 
number of independent variables to 3; in this case, they 
are chosen as n^, . The entropic terms due 

to a loss of translational entropy and due to combinato- 
rial gains by three-phase coexistence are constructed in 
complete analogy to eqs. (46) and (47). 



C. Fractionated three-phase equilibrium conditions 

We minimize the fractionation free energy presented 
in the last subsections with respect to the volume frac- 
tion and sequence distribution of the emerging shadow 
phase (s). Insertion of the free-energy densities of the dif- 
ferent states with fractionation into eq. ( 48 1 and subse- 



quent differentiation of /f rac with respect to the variables 

(21 (2) (2) fll (1) (1) • 

v y ' , n 2 , rin or v K > , n 2 , n\ give equation systems 



9/frac <9/frac 



P =1.2 



(53) 



exemplified in appendix|Dj Solutions are obtained numer- 
ically with a Newton-type procedure (cf. appendix |E[) . 



D. Three-phase transition lines 

Upon gradually decreasing or increasing the incompat- 
ibility L\ at fixed A in the three-phase state, boundaries 
of the three-phase region, x at A > A c , and \^ (at 
A < A c with our simplified lamellar free energy), are in- 
dicated by a zero of the free energy eq. ( [48] ) due to frac- 
tionation: Either the minority phase's volume fraction 
tends to zero (characteristic of a shadow), its sequence 
concentrations approach the A-defined ones, or its order- 
parameter amplitude tends to zero. Analysis of eq. ( 53 ) 



shows that in our system the first alternative is realized, 
which simplifies the set of equations for the transition 
lines. 



IV A 



In the case of sec 
tropic contributions eqs. (146b , (|47j) to /f r , 
fraction v^ 2 ' of the lamellar shadow phase yields 



an expansion of the en- 
in the volume 



A/ V ol. red. + A/, 



comb. 
,(2) 



(54) 



,(2) 



1 3 / T7 (2) 




o 



Similarly, one can expand the deviations from A-defined 
concentrations in the two-phase cloud state: 

-p v {\) = (p v {\)-n^)v^ (55) 

+(ma)-4 2 >) (v^y+o^y^ 

Hence the lowest-order term of /f rac is linear in v^ 2 \ 

A/ frac (^ 2 \r4 2) ,r4 2) ) (56) 
= a(n 2 2 \nf ) )v^+o((v^f 



and the coefficient a must be minimized in order to de- 
termine x^ 1 "* and the sequence distribution of the shadow 
phase. The phase transition line from the global lamellar 
to the fractionated three-phase state can be treated in 
complete analogy by taking the limit yW — > 0. 



V. PHASE DIAGRAMS 

In the following, we present the lines of macroscopic 
and lamellar phase separation and the boundary lines 
of three-phase coexistence obtained from the minimiza- 
tion of the fractionation free energy. The critical line of 
macroscopic phase separations of the disordered melt is 
well known already from approaches based on the multi- 
component picture [7J [S]- Also the discussion of the 
pure microphase separation transition within mean-field 
theory can be found elsewhere [3l EU [29]. Our focus 
here is on the coexistence of homogeneous and lamel- 
lar phases with fractionated sequence distributions. The 
point (A c , Lxc) where the transition curves from the dis- 
ordered toward macroscopic, respectively lamellar, phase 
separation meet will be mostly referred to as a multi- 
critical point without further classification. Partitioning 
of sequences will be shown via distribution diagrams in 
sec. |VII[ in comparison with SCFT calculations. 



A. Triblocks with small M 

To exemplify the phase behavior of triblocks with 
M < 7 segments per block, we discuss the results for 
M = 3 shown in fig. [3] To explore the emergence and 
growth of the various phases, we follow the path indi- 
cated by arrows in the plot, starting at a block corre- 
lation A > A c = — !/3: The first instability of the dis- 
ordered melt is toward homogeneous phase separation, 
indicated by the peak at zero wave number of the global, 
second-order structure factor (cf. the solid curve in the 
bottom inset). Upon increasing incompatibility L\ (bot- 
tom vertical arrow), the dotted line (x ) marks the on- 
set of three-phase coexistence via a fractionated lamel- 
lar shadow phase with volume fraction = 0. (A 
fractionated lamellar shadow was already predicted by 
Monte Carlo simulations [TH].) This fractionated phase 
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FIG. 3. (Color online) Phase diagram for triblock copolymers 
with M = 3. Solid (red) line marks macroscopic (two ho- 
mogeneous, A- and B-rich phases), dashed (blue) line marks 
lamellar phase separation (order-disorder-transition ODT) of 
the disordered state. Dotted (orange) line: onset of three- 
phase coexistence, at which the two homogeneous phases are 
the cloud and a fractionated lamellar phase shadow emerges; 
dot-dashed (green) line is the lamellar cloud boundary. A 
circle marks the multicritical point. Bottom inset: second- 
order structure factor in the global A-distribution, at A = 
(solid), at the critical correlation A c = — V 3 (dotted), and at 
A = —0.5 (dashed). In this and the following plots of this 

part, the length scale is ((n ■ R,,) 2 ^ 1 ^ 2 — R,,/\fd. Top inset: 
sketch of three-phase coexistence. 



sets in with finite amplitude, and with finite wave num- 
ber, since its copolymer-enriched sequence distribution 



(see fig. 16 below) causes the structure factor to be dif- 
ferent from the global one. On further increase of the 
incompatibility (along the top vertical arrow), the lamel- 
lar volume fraction grows. Now, keeping L\ constant, 
and proceeding toward smaller values of A (following the 
horizontal arrow), the volume fraction of the lamellae in- 
creases further. Finally, at some A < A c , one reaches 
the end of the three-phase coexistence (indicated by the 
dot-dashed line), and lamellae take over to be the cloud 
phase with volume fraction i/ 2 * 1 = 1. Consistently, start- 
ing at A < A c and small incompatibilities, the disordered 
melt undergoes lamellar phase separation (at the incom- 
patibilities on the dashed line) due to the peak of the 
A-defined structure factor at a finite wave number (cf. 



eq. (391 and the dashed curve in the bottom inset). With 



our simplified free energy eq. ( 26 ) , via which the insta- 



bility toward a global lamellar phase rests solely on the 
fc-dependence of this second-order structure factor, the 
lamellar cloud boundary of three-phase coexistence is al- 
ways located in the half-plane A < A c . Upon crossing 
the dot-dashed boundary line from this side, two addi- 
tional homogeneous phases with homopolymer-enriched 
sequence distributions appear as shadows. 

As visible in fig.[3| three-phase coexistence prevails in a 
large parameter region. However, since our lamellar free 



energy is limited to small order-parameter amplitudes, 
the results may be unreliable at very large values of the 
incompatibility. An alternative scenario would be global 
lamellar phase separation at higher L\ (see fig. 14 below) . 

At the critical correlation A c , the maximum at fcrj = 
of the global structure factor broadens (see dotted curve 
in the bottom inset in fig. [3| , announcing the continuous 
growth of the optimal wave number from zero when low- 
ering A. Qualitatively, we observe this transition from 
global macroscopic to global lamellar phase separation 
for all random triblock copolymers with M < 7 (cf. the 
case discussed before eq. (30)), while the exact position 
of the Lifshitz point (A c , L% c ) depends on M. This point 
of diverging lamellar wavelength limits the three-phase 
region toward low incompatibilities. 

The lamellar wave numbers on the boundary lines of 
fractionated three-phase coexistence as a function of A 
are displayed in fig. |4j Note that the simplified free en- 
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FIG. 4. (Color online) Lamellar wave numbers for triblocks 
with M = 3. Dot-dashed (blue) line: global lamellar phase at 
A < A c , between the order-disorder transition and the onset 
of three-phase coexistence with v^ 2 ' = 1 (hatched region of 
fig. M). Dotted (orange) line: fractionated lamellar shadow 

„(2r_ 



0). Inset: enlarged around the multicritical point. 



ergy for microphases (see eq. (39)), predicts that at a 
given A, the wave number of global lamellar phase sep- 
aration (hatched region in fig. [3| does not change with 
increasing L\. The lamellar wave number can be shifted 
only due to fractionation, i.e., by an increased content 
of alternating sequences. We find that, on increasing L\ 
in the three-phase region, the fractionation and thereby 
the wave number in the lamellae increase, i.e., the lamel- 
lar spacing decreases. This is in agreement with findings 
for global microphase separation in random copolymers 
within mean-field theory [2J. 

The wave number of fractionated lamellae vanishes at 
the Lifshitz point (A c ,Xc), as does the wave number of 
global lamellar phase separation. The inset in fig. [4] 
shows the behavior of the fractionated wave number in 
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the vicinity of the Lifshitz point. For A < A c , the three- 
phase region can be entered at two different incompatibil- 
ities, with different wave numbers of the lamellar shadow. 
A closer look is cast onto this remarkable feature of the 
phase diagram in the detail of the boundary lines and a 
map of the lamellar phase's volume fraction around the 
Lifshitz point in fig. [5] The line of fractionated lamel- 
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FIG. 5. (Color online) Volume fraction of lamellar phase and 
three-phase coexistence lines around the multicritical point 
(A c , L\ c ) for triblocks with M = 3. Line styles as in fig. [3] 

lar shadows displays a reentrant behavior, especially it 
does not reach the Lifshitz point for A \ A c , but via a 
spiraling path invading the region A < A c . Except for a 
very small region of the parameter space, fractionation 
suppresses global lamellae with diverging wavelength in 
the vicinity of the Lifshitz point, in favor of, first, macro- 
scopic phases and, at higher incompatibilities, fraction- 
ated lamellae with finite wavelength. 
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FIG. 6. (Color online) Lamellar order-parameter amplitude 
along the boundaries of three-phase coexistence for M = 3. 

The scaling of the order-parameter amplitude on ap- 
proach to the multicritical point along the transition lines 



to three-phase coexistence (A /* A c ) is shown in fig. 6] 
The amplitudes of fractionated lamellar shadows (on the 
dotted line in fig. [5j in the range A < A c ) are marked by 
open diamonds, those of global lamellar (cloud) phases 
(on the the dot-dashed line in fig. [5]) by solid triangles, 
those of the coexisting homogeneous shadows by open 
squares. According to the fit performed to the latter 
case, the amplitudes vanish linearly in AA := |A — A c | 
at the Lifshitz point (the same exponent is found for the 
lamellar cloud amplitude). 

In order to analytically extract the exponent of the 
order-parameter amplitude in the vicinity of the multi- 
critical point, we solve t he eq uation of the lamellar cloud 
line i/ 1 ) = (see sec IV D I for the deviations of the 
sequence concentrations in the fractionated macroscopic 
shadow phases from the global ones, An„ := riu —p„(\), 
with a power series ansatz 



An l ,(AA) = ^c^(AA) J 



(57) 



The series' coefficients of the equation system in AA 
can be calculated for Me [3, . . . , 6], cases in which the 

1/2 

wave number ko of global lamellae vanishes oc (AA) at 
(A c , Xc) (as expected for a Lifshitz point [3])- For M = 3, 
consistent expansion up to (AA) 4 yields, along the lamel- 
lar cloud boundary line, 

An, = - 144 ^(AA) 2 +G((AA) 3 ), v = 2,3. (58) 
55 

When inserting these dependencies into expansions of 
the optimal wave number, the structure factor etc. (cf. 
eq. (p8|)), we indeed find the critical exponent 1 for the 
amplitude <j m in the lamellar cloud phase, 



oc AA, A /■ A c . 



(59) 



Moreover, the slopes of the transition lines I<x(A) from 
the disordered to the global lamellar state (xm(A)) and 
from global lamellae to fractionated three-phase coexis- 
tence (x' 2 '(A)) can be shown to be equal at (A c ,Lxc) : 



X (2) (A)-Xm(A)cx(AA) 2 , \ /• \ c 



(60) 



Continuous triblocks 



Representative of triblocks with M > 7 segments per 
block, the phase diagram for continuous random triblock 
melts is shown in fig.[7| Again, for A > A c , the dotted line 
marks the emergence of a lamellar shadow in addition 
to the two homogeneous phases. The lamellar volume 
fraction grows with increasing L\ and with decreasing 
A. On the dot-dashed line, the lamellar phase takes over 
to be the cloud phase and coexists with two fractionated 
homogeneous shadows. 

In comparison to the case M < 7 (see fig. [3]), the 
three-phase coexistence region seems to be larger. (Still, 
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FIG. 7. (Color online) Phase diagram for continuous tri- 
blocks. Line styles as in fig. [3] Crosses indicate the end 
points of the lines of metastable, global phase separations, 
macroscopic for A < A c (x) and lamellar for A > A c (+). Bot- 
tom inset: global second-order structure factor at A = —0.45 
(solid), at the critical correlation A c = —0.464 00 (dotted), 
and at A = —0.47 (dashed). Length scale as in fig. [3] 



the predictions are restricted to incompatibilities that do 
not exceed considerably those of the order-disorder tran- 
sition.) The multicritical point is not only located at 
a smaller critical block correlation A c and a higher in- 
compatibility, but is also qualitatively different: As dis- 
cussed in sec. |IIIB| the wave number of the first global, 
ordered structure (when starting at low incompatibili- 
ties in the disordered state) is discontinuous at A c for 
M > 7. Thus when reaching A c from above, the morphol- 
ogy of the ordered phase changes from two homogeneous 
phases (zero wave number fcn = 0) to one lamellar phase 
with finite wave number fcrj,c- This feature is revealed 
in more detail in the plot of lamellar wave numbers in 
fig. [8] At the multicritical point, the lamellar wave num- 
bers in the fractionated state also tend to the finite value 
ko.c = 0.3267T. The wave number in the fractionated 
lamellar shadow attains a slightly smaller, minimal value 
at a correlation A > A c (cf. the top inset in fig. [8]). Due 
to the two peaks of the global structure factor around 
multicriticality (see the inset in fig. [8b, metastable global 
lamellae occur in a small range of block correlations, 
—0.461 23 > A > A c , where the free-energy functional's 
absolute minimum indicates global macroscopic phase 
separation. Inversely, global macroscopic phase separa- 
tion persists as a metastable state for —0.5 < A < A c (at 
A = —0.5, the curvature of s(k 2 ) at k = changes accord- 
ing to eq. (|30|). These metastable transition lines, whose 
end points are hardly resolvable in fig. [7J are displayed 
in fig. [9j together with the actual transition lines and a 
map of the lamellar volume fraction around the multi- 
critical point (note the zoom to an even smaller region 
than in fig. |5]). On increasing incompatibility from the 
fractionation onset, the lamellar volume fraction grows 
(for A > A c ) or decreases (for A < A c ) rapidly to level out 
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FIG. 8. (Color online) Lamellar wave number for continu- 
ous triblocks. Dot-dashed (blue) line: global lamellar (cloud) 
phase at A < A c ; dotted (orange) line: fractionated lamellar 
shadow (?/ 2 ' = 0); solid (cyan) line: metastable global lamel- 
lae due to a second peak of the structure factor s(fe 2 ), shown 
at A = —0.4625 (triangle) in the bottom inset. Circles mark 
the wave numbers of the coexisting states at the multicriti- 
cal point, crosses mark the end points of metastable global 
lamellar/macroscopic phase separation lines. 



11.48 




E 1 1 .42 

o 

o 



11.38 - (K, L Xc ) = (-0.4640, 1 1 .431 6) 



-0.466 -0.464 -0.462 
block correlation X 



-0.46 



FIG. 9. (Color online) Volume fraction of the lamellar phase 
around the multicritical point (A c , L\c) for continuous tri- 
blocks. Boundary lines of three-phase coexistence as in fig. [7] 
Additional thin lines mark metastable, global phase separa- 
tions: solid (cyan) for A > A c , with end point marked by a 
cross: lamellar phase; dashed (red) for A < A c : homogeneous 
phases. 



at a value of about 0.6. At multicriticality, the transition 
lines from the disordered to the global lamellar state and 
from global lamellae to fractionated three-phase coexis- 
tence differ in their slopes, in contrast to the behavior at 
the Lifshitz point. 

Despite the discontinuity of the wave number at the 
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FIG. 10. (Color online) Scaling of the lamellar order- 
parameter amplitude along three-phase boundaries for con- 
tinuous triblocks. 



critical correlation A c for M > 7, the boundary lines 
of fractionated three-phase coexistence are single-valued 
around the multicritical point. Hence, in this case, we 
can determine numerically the critical exponent for the 
decay of the lamellar order-parameter amplitude along 
both boundary lines (see fig. 10). The exponent 0.5, 
found along both lines, is reminiscent of mean-field be- 
havior. Note that for triblocks with M < 7, we derived 



a different critical exponent, viz. 1 (cf. eq. (59)) 



C. Fractionation with restored wave number 
dependence 

In this section, we aim at testing the fractionation 
scenario with the complete fourth-order expansion of 
the Landau-type free energy for structured phases, in- 
stead of the simplified version eq. (26). To this end, 



we accounted for the wave number dependence of the 



fourth-order terms of the functional eq. ( 23 ) in eq. ( 40 ) 
in sec. MI C 21 The effects of the wave number varia- 
tion within our fractionation scheme can be observed in 
fig. |11[ for random continuous triblocks. The bound- 
ary between global macroscopic phase separation and the 
three-phase region at A > A c is located at lower incom- 
patibilities than that obtained with the simplified free 
energy (cf. fig. [7]). Global lamellar phase separation is 
found to be stable in a larger parameter region and to 
extend into the half-plane A > A c . However, upon fur- 
ther increasing x i n the system with global lamellar phase 
separation at A > A c , we find a reentrance into the frac- 
tionated three-phase coexistence. Note that the ampli- 
tude of the lamellar shadow at the onset of fractionation 
attains a reasonably small value also at a block correla- 
tion distant from the critical one (cf. the sinusoidal profile 
in fig. 111]). 



The main advantage of the lamellar free energy eq. ( 40 ) 



is the principal possibility of global lamellae also at A > 
A c , since the optimal wave number changes with increas- 
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FIG. 11. (Color online) Phase diagram for continuous tri- 
blocks with ^-dependence of the free-energy functional eq. 
(23 1. Line styles as in fig. [3] Crosses indicate the end points 
of the line of metastable macroscopic phase separation for 
A < A c (x) and of metastable lamellae for A > A c (+). Top 
right inset: profile of A fraction in lamellar shadow at A = 0. 



ing x even in a fixed sequence distribution (similar to 
the mechanism of global microphase separation invoked 
in ref. |24j . which, however, considered one-component 
diblock copolymers only). Global lamellar phase separa- 
tion is found to follow the three-phase coexistence at high 
incompatibilities also within SCFT (see fig. 14 below). 

The scaling of the order-parameter amplitude on ap- 
proach to the multicritical point is exctracted from the 
regularly shaped lamellar shadow line in fig. [12] Both 
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FIG. 12. (Color online) Order-parameter amplitude of ho- 
mogeneous clouds and lamellar shadow at the onset of frac- 
tionation for continuous triblocks, with fc-dependent fourth 
order. 



the lamellar shadow and the macroscopic cloud ampli- 
tudes vanish with an exponent of 0.5, corroborating the 
findings with the simplified lamellar free energy. 
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VI. ALTERNATIVE APPROACH: SCFT 
A. Method 

An alternative method to determine the phase behav- 
ior of random triblocks employs self-consistent field the- 
ory (SCFT) [301 OH]- In order to analyze the phase co- 
existence of homogeneous and lamellar phases with finite 
volume fractions, it starts out from the grand-canonical 
partition function, 



^=£n 



1 / QoC 



{N v } v 



N v \ \LQ Q 



(61) 



J V [r{s)] exp {-H x -H K - H w } , 

with In Q v being the excess chemical potential of species 
v, v = 1, . . . , 6. Q denotes the configurational partition 
function (without translation) of a single, non-interacting 
Gaussian chain. Via the incompressibility demand (see 
below), the sum over the sets {N v } of species numbers 
is restricted by the constraint, $^„ =1 N v = N . There- 
fore not all the fugacities are independent, and we set 
C,aba = 1- By virtue of the symmetry A^B of the A- 
distribution and the coexisting phases, C,aaa — Cbbb, 
(aab — Cbba, and (aba — Cbab- 

Similarly to the formalism in the previous sections, A- 
and i?-density fields with their respective auxiliary fields 
wa and Wb are introduced to decouple the interacting 
chains. The incompressibility constraint is accounted for 
by an additional Lagrange field £ and automatically im- 
poses the constraint on the species numbers. Within the 
saddle point approximation, we obtain the excess grand- 
canonical potential, g = G/N + 1, per molecule: 

g = ^ J d 3 r {xL(f>A4>B ~ w a 4>a - w b 4*b} , (62) 

where the saddle point values of the fields and densities 
are determined by the self-consistent set of equations 



" <f>B = 1, 

w A = xL<Pb + f , 
wb = xL<Pa + £, 



i(r) = -j2CuV 



5Q, 



^-f 8w B (r) 
The global concentration p v of species v is given by 

Pu = ^ J d 3 r (f> v (r) = CuQu- 



(63a) 
(63b) 
(63c) 



(63d) 
(63e) 



(64) 



The saddle point equations involve the partition func- 
tions, Q v , of single copolymer chains of species v in the 
external fields wa and wb- 



exp 



ds w v (r(s), s) 



with 



(65) 



w v (r(s), s) := 



1 



w A {r(s)) + w B (r(s)) 



with the conformational average defined in eq. ( |15[ ). 
In the following, only the continuum limit of Gaussian 
chains is considered. For a structured phase, the Q v 
and density profiles are expressed in terms of statisti- 
cal weight propagators q„(r, s), ql(r, s) along a Gaussian 
chain, 



q v {r,s) = (exp<^ - / ds' w v (r(s'), s') 



Ql{r,s) 



exp 



ds' w v (r(s'), s) 




(66a) 



(66b) 



where and H^~ r % )=r) 



W(t-(s) = t-) 

are the conformation 
statistical weights for a chain of length s having its end 
point at r and for a chain of length L — s having its 
start point at r, respectively. The single-chain partition 
functions Q v are calculated according to: 



Q v = 



d d rq u (r,s)qUr,s), Vsg[0,L]. (67) 



0, 



The propagators obey the modified diffusion equations: 
d_ 
ds 
d_ 
ds 



A r + w„j q„(r,s) 
A r - w u J ql(r,s) = 



(68a) 



(68b) 



These partial differential equations are solved via a spec- 
tral method [32] ■ As a result, we obtain the equilibrium 
spacing and the free energy of the lamellar phase, as well 
as detailed composition (concentration) profiles of the 
different species in a lamellar domain. An example of a 
composition profile is shown in Fig. [13] for A = at the 
lamellar cloud point, L\ — 9.389 19. 

The canonical free energy can be obtained via a Leg- 
endre transformation: 



F = G + N v In C, v + N In L '" 



LQ 



(69) 



Thus the excess Helmholtz free energy / per molecule 
takes the form 



f = 



F 

N 



In 



Go 



Vv {\k Vv - 1) + ^ J d 3 



r <Pa<Pb 



(70) 



Pu ln Q v 



V 



d 3 r {wA't'A + w B (^b) 
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FIG. 13. (Color online) Local composition (A-segment den- 
sity) profiles of the lamellar cloud phase at A = 0. Solid line: 
total density, symbols (key in the plot): due to one copoly- 
mer species. The spatial coordinate x is normalized by the 
lamellar spacing, A = 4.148i? c , where R c denotes the rms 
end-to-end distance of a non-interacting triblock copolymer. 



The first term corresponds to the entropy of mixing of 
the different species, the second term quantifies the free 
energy due to the repulsion of unlike monomer types, and 
the last two terms describe the loss of conformational 
entropy of the polymers in a spatially inhomogeneous 
environment. 
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FIG. 14. (Color online) Phase diagram for continuous tri- 
blocks within SCFT. Line styles as in fig. [3] Symbols highlight 
data: triangles (green) lamellar cloud points, squares (blue) 
order-disorder transitions. The inset shows a detail, with the 
three-phase boundaries from the analytical method (A) added 
(cf. fig. Ill, in dashed (orange) homogeneous, in dot-dashed 
(green) lamellar cloud points. Thin dots (green) extrapolate 
the SCFT lamellar cloud points to the multicritical point. 



—0.17 < A < 0.43, the data were calculated with a spa- 
tial resolution of 32 Fourier components, in the remaining 
range with 12 components. 



B. Three-phase coexistence lines and fractionation 

1. Homogeneous cloud phases 

If we approach three-phase coexistence by increasing 
the incompatibility L\ from a low value at fixed A, 
the lamellar phase (shadow) will emerge from the ho- 
mogeneous, A-rich and B-rich phases (clouds) with an 
infinitesimal volume fraction. At the onset of three- 
phase coexistence, the sequence distribution of the two 
homogeneous cloud phases is a A-defmed one. In the 
grand-canonical ensemble, we determine the two inde- 
pendent excess chemical potentials, Caaa and (aab, of 
the cloud phases as to reproduce the composition of 
the A-distribution. Since the incipient lamellar phase 
can exchange polymers with the cloud phases, its prop- 
erties are calculated in the grand-canonical ensemble. 
To this end, we minimize the grand-canonical poten- 
tial, g, at given Caaa and (aab with respect to the 
lamellar period or spacing A. The onset of three-phase 
coexistence occurs at the incompatibility, at which the 
so-minimized grand-canonical potential of the lamellae 
equals the grand-canonical potential of the cloud phases. 
The (A, x)-points, at which the homogeneous, A-rich and 
B-rich phases are the cloud phases, are shown as a dot- 



ted curve in the phase diagram fig. 14 In the range 



2. Lamellar cloud phase 

As we progress into the three-phase coexistence to- 
ward larger incompatibilities L\, the volume fraction 
of the lamellar phase grows, while that of the homoge- 
neous, A-rich and i?-rich phases decreases. At the end 
of three-phase coexistence, the lamellar phase occupies 
the entire volume, and the homogeneous phases contin- 
uously disappear with a vanishing volume fraction. In 
order to determine this cloud point of the lamellar phase, 
we calculate the properties of the latter in the canonical 
ensemble, where its sequence distribution is fixed to A- 
distribution. The canonical free energy / is minimized 
with respect to the lamellar spacing A. Then, the two 
independent excess chemical potentials for this optimal 
lamellar structure are measured, and the properties of 
the incipient homogeneous phases are calculated in the 
grand-canonical ensemble at the so-determined chemical 
potentials. Finally, L\ is adjusted such that the lamel- 
lar cloud and the incipient homogeneous shadow phases 
have the same grand-canonical potential at identical ex- 
cess chemical potentials. The resulting boundary points 
of three-phase coexistence toward large L% are marked in 
fig. [14] by triangles on a dot-dashed line. Twelve Fourier 
components were considered in this calculation. 
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C. Phase coexistence with finite volume fractions 

The properties of a general fractionated state of three 
coexisting phases are computed in the grand-canonical 
ensemble. As in sec. VI B 2| we consider a lamellar phase 
(marked by the superscript (2)) with volume fraction 
v^ 2 \ coexisting with two homogeneous, A-rich and i?-rich 
phases, with joint volume fraction 1 — v^ 2 \ A fractionated 
state with given volume fractions is located by simulta- 
neously adjusting the two independent excess chemical 
potentials Caaa an d Caab and the incompatibility L\ 
such that the weighted sum of the sequence concentra- 
tion p^ , respectively p^ , in the lamellar phase and p^ , 

respectively p^\ in the two homogeneous phases gives 
the global concentration of a A-distribution (cf. eqs. ( 29 ) 



and (50)), and such that the grand-canonical potentials 
of all three phases are equal [33] . — g( 2 \ In the limit 
w ( 2 ) _^ o, we recover the cloud points of the homogeneous 
phases, in the limit i/ 2 ) -> 1 we recover the cloud points 
of the lamellar phase. In contrast to the phases at their 
cloud points, none of the coexisting phases with finite 
volume fractions displays a A-distribution (cf. sec 
below). 
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FIG. 15. (Color online) Incompatibility L\ against volume 
fraction v' 2 ^ of the fractionated lamellar phase at A = 0. The 
end points of the curve mark the limits of three-phase coex- 
istence, at which the lamellar phase is the incipient shadow, 
i/ 2 ' = 0, respectively the cloud, t;' 2 - 1 = 1. The top inset shows 
the lamellar A density profiles in these limits. The bottom 
inset shows the spacing A of lamellae in the three-phase coex- 
istence region and that of global lamellae for L\ > 9.389 19. 



The gradual change of the volume fraction of the lamel- 
lar phase upon increasing the incompatibility L\ at A = 
is shown in fig. |15| The inset presents the composition 
(A-segment density) profiles of the lamellar phase at its 
shadow point L\ — 8.820 43 (dashed line), and at its 
cloud point L\ = 9.389 19 (solid line). We observe that 
the lamellar shadow's profile, though it is not confined to 
a single harmonic, matches quite well the profile obtained 
from the analytical method (see the inset in fig. 11), es- 



pecially in the amplitude. In contrast to results for one- 
component diblock copolymer melts [24 , but in agree- 
ment with predictions of random phase approximation, 
the lamellar spacing decreases upon increasing L\- 



VII. FRACTIONATED SEQUENCE 
DISTRIBUTIONS 

In this section, we invoke both the analytical and the 
SCFT method to obtain detailed sequence distributions, 
which show the fractionation or sequence partitioning ac- 
cording to the coexisting phases' structures in random 
continuous triblocks. In figs. [16] and [T7j the sequence 
distributions of the coexisting phases are presented by 
means of composition triangles: Each corner represents 
one of the sequence classes defined in eq. (28), a point 
within the triangle one sequence distribution. Due to 
the A^B exchange symmetry of species combined into 
one sequence class, the distributions of the two homoge- 
neous phases within a macroscopically phase-separated 
state coincide in this triangle. 
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FIG. 16. (Color online) Sequence distribution triangle for ran- 
dom continuous triblocks at various block correlations, with 
the extended analytical method, cf. the phase diagram in 
fig. 11 The diagram's center (+) corresponds to equal concen- 
trations of all sequences (p„ = 1/3, v = 1,2, 3). One sequence 
distribution is represented as a linear combination of the vec- 
tors pointing from the center to the corners, each vector scaled 
with the concentration deviation 3 (p v — 1 /3) /2. Distributions 
defined by A lie on the (red) curve, with A ranging from — 1 at 
the triangle's bottom left corner to +1 at its top. Solid sym- 
bols on this curve mark the sequence distribution of the cloud 
phase(s) at the boundary line(s) of three-phase coexistence. 
Off-curve solid symbols mark the distributions of the coexist- 
ing shadow phases. Open symbols display the distributions of 
the coexisting states at equal volume fractions (i/ 2 - 1 = 0.5). 

In fig. |16[ we present the fractionated distributions 
obtained by the analytical method with the restored k- 
dependence of the lamellar free energy. The sets for three 
supercritical values of the block correlation, A > A c , A = 
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0.2 (diamonds), A = (circles), and A = —0.2 (up tri- 
angles), visualize the following fractionation mechanism: 
On the curve of A-distributions, the solid symbol indi- 
cates the sequence distribution of the homogeneous cloud 
phase(s) at the onset of fractionated three-phase coexis- 
tence. The solid symbol of the same shape and color to 
the bottom right of the curve, marks the distribution of 
the coexisting lamellar shadow phase (with zero volume 
fraction). The finite deviation of the lamellar shadow's 
sequence distribution from the A-distribution shows that 
the transition to three-phase coexistence is discontinuous. 
Upon increasing incompatibility, the lamellar phase's vol- 
ume fraction increases (cf. fig. [9]), and its sequence dis- 
tribution departs ever more from the A-distribution (the 
open symbols to the bottom right of the A-curve dis- 
play lamellae at 0.5 volume fraction). Sequence class 2 
(AAB/BBA) substantially accumulates in the lamellar 
phase, also class 3 (ABA/ BAB). Moreover, since the ra- 
tio of these two sequence concentrations differs from the 
A-defined ratio P2W/P3W, the fractionated sequence 
distribution in the lamellar phase does not ensue from 
merely expelling homopolymers into the coexisting ho- 
mogeneous phases at a constant ratio of the other two 
sequence classes. As the volume fraction of the homoge- 
neous, initial cloud phase(s) decreases, their distribution 
(at volume fraction 0.5 marked by open symbols to the 
top left) deviates increasingly from the A-curve, showing 
in turn a particular depletion in AAB/BBA sequences. 

For A = —0.4, the reentrant behavior of the three- 
phase boundary line, cf. 
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BBB 



fig. 11 



gives rise to various co- 
existing distributions (down triangles). Upon increas- 
ing x m the two homogeneous phases, the three-phase 
region appears with a lamellar shadow (nearly on the 
curve of A-distributions, shift to the bottom hardly visi- 
ble), which grows with x i n volume fraction until it be- 
comes a lamellar cloud (now the symbol on the curve of 
A-distributions) coexisting with two homogeneous shad- 
ows (triangle shifted slightly to the top). The homoge- 
neous shadows are homopolymer-enriched and depleted 
in alternating sequences. At an even higher x> the global 
lamellar phase gives way to a three-phase coexistence 
again. The topmost triangle represents the distribution 
of the homogeneous shadows at this reentrance. For 
A < A c , this lamellar cloud line is the only three-phase 
boundary. The topmost symbols for the critical and sub- 
critical correlations A = A c and A = —0.5 show the dis- 
tributions of the coexisting homogeneous shadows which 
deviate markedly from the A-distributions. 

In the distribution triangle of fig. [TTJ we present the 
SCFT results for the sequence distributions of the co- 
existing phases for A = 0, 0.25, and 0.5. Again, the 
distributions of the cloud phases are represented by solid 
symbols on the solid curve of A-distributions. For each 
value of A, the distributions at the beginning and the 
end of three-phase coexistence are shown. At the lower 
incompatibility, the homogeneous phases are the clouds, 
and the coexisting lamellar shadow corresponds to the 
respective solid symbol shifted to the lower right corner, 
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FIG. 17. (Color online) Sequence distribution triangle for 
random continuous triblocks within SCFT, with the distribu- 
tions of the coexisting phases at the beginning and at the end 
of three-phase coexistence for A = (blue circles), 0.25 (ma- 
genta squares), and 0.5 (red stars), respectively. For A = 0, 
open circles mark the distributions of the coexisting phases 
at Lx = 9.017 23, where the lamellar phase comprises half of 
the volume. The solid line represents A-distributions. 



with its distribution enriched in AAB/BBA sequences. 
At the higher incompatibility, the lamellar phase occu- 



pies the total volume, v 



(2) 



1, and its distribution is 



represented by the cloud symbol on the A-curve. The 
distribution of the coexisting homogeneous shadows cor- 
responds to the symbol shifted to the upper left side of 
the triangle. Two open circles mark the distributions 
for equal volume fractions of the macroscopic and the 
lamellar phase-separated state, — 0.5, at A = 0. In 
this situation, none of the coexisting phases is character- 
ized by a A-distribution; the homogeneous phases are rich 
in homopolymers, while alternating sequences segregate 
into the lamellar phase. In comparison to the analytical 
results for the distributions at A = 0, apart from the qual- 
itatively different feature of a lamellar cloud at higher in- 
compatibilities, the sequence fractionation is found to be 
weaker. Note, however, the smaller transition incompat- 
ibilities to three-phase coexistence within SCFT, which 
also result in smaller order-parameter amplitudes. 

Detailed sequence distribution diagrams for lamellar 



and macroscopic phases at A = are displayed in fig. 18 
for the analytical method, and in fig. [19] for SCFT. The 
representation of all six species' concentrations addition- 
ally visualizes the segregation within a sequence class 
into A- and £?-rich subspecies between the two homoge- 
neous phases, which allows for an estimate of the macro- 
scopic A excess amplitude at different stages of fraction- 
ated three-phase coexistence. Due to the homogeneous 
phases' A ^ B exchange symmetry, only the distribu- 
tion of the v4-rich, homogeneous phase is shown. (The 
chart for the B-rich phase looks the same as the one de- 
picted for the A-rich phase, only with letters A and B ex- 
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3 BBA tzszsz BAB 
AAA gzzza AAB kxxxe ABA 




X-distrib 
(X = 0) 



macr. J$ m -> 
A-rich v ( ' = 0.5 



FIG. 18. (Color online) Detailed sequence distributions of the 
coexisting phases at A = from analytical method (sec. V C I, 
cf. fig. |16| Leftmost chart: A-distribution of the disordered 
melt. Pairs of charts: distributions of the A-rich, homoge- 
neous and the lamellar phase; left: at the onset of three-phase 
coexistence; right: at a lamellar volume fraction of 0.5. 



BBB i/ / / BBA i vvn BAB 
AAA f///// AAB r vww i ABA 




Lx = 8.82 




Lx = 9-39 




X-distrib. 
(X=0) 



macr. 
A-rich 



lam., rr 
v® = A 



macr. lam., 
rich v ( ' = 1 



FIG. 19. (Color online) Detailed sequence distributions of 
the coexisting phases at A = obtained with SCFT. Pairs 
of charts: distributions of the A-rich, homogeneous and the 
lamellar phase; left: at the onset; right: at the end of three- 
phase coexistence. 



changed in the key.) The distributions obtained by both 
methods agree well. While both diagrams reveal the pref- 
erence of the fractionated lamellar phase for AAB / BBA 
sequences, the accumulation is more distinctive in the an- 
alytical results, already at the onset of fractionation (cf. 
the central charts). Corresponding to the higher onset 
incompatibility, the macroscopic segregation into A- and 
£?-rich subspecies also is at a more advanced stage. 



able to analyze accurately and in detail the vicinity of the 
multicritical point (A c ,Xc), whose quality is found to de- 
pend sensitively on the number M of segments per block. 
For a small number of segments per block (M < 7), the 
wave number of the global ordered state grows contin- 
uously from zero, when decreasing A from A c . A reen- 
trance into the fractionated three-phase coexistence is 
observed for A < A c . The critical exponent for the lamel- 
lar order-parameter amplitudes on approach to (A c , Xc) 
is 1, along both three-phase coexistence boundaries. For 
more segments per block (M > 7), the structure factor 
of the A-distribution develops a second peak at finite k, 
such that the wave number of the global ordered state is 
discontinuous at A c . Contrasting with the case M < 7, 
we find a critical exponent of 0.5 for the amplitudes of 
both lamellar and homogeneous phases, along both three- 
phase boundaries. This behavior might be due to the 
intersection of transition lines to metastable, global or- 
dered phases at multicriticality. 

With the simplest version of the free-energy functional, 
cq. (26), a (global) lamellar cloud phase with A-dcfincd. 
concentrations can occur only at A < A c (for our sys- 
tem, a result qualitatively different from the SCFT pre- 
dictions; see sec. | VIIIB| below). An enhanced version 



of our theory abandons this restriction by restoring the 
wave-number dependence of the quartic vertices in the 
lamellar free-energy function, eq. ( 40 1 ; see the location 



of the lamellar cloud boundary for random continuous 
triblock copolymers in fig. |11| The critical exponent 
of 0.5, found for the order-parameter amplitudes along 
three-phase coexistence lines with the simplified theory, 
is corroborated by the enhanced analytical theory. 



With the complete wave- vector dependence of eq. ( 23 1 , 



at fixed A, a global lamellar phase can attain a lower free 
energy than global macroscopic phase separation at an 
incompatibility Xh > Xh ~ a mechanism of microphase 
separation proposed by Leibler and co-workers O 124) . 
Via our parameterization of a fractionated three-phase 
coexistence, we take into account more degrees of free- 
dom and find, instead of this mechanism, a refined com- 
petition to be effective: A structured phase first becomes 
stable in a subsystem with vanishing volume fraction and 
with a sequence distribution different from the global one. 
This onset of three-phase coexistence indeed occurs at a 
smaller incompatibility x < Xl than that of the global 
microphase separation conjectured by Leibler. 



B. Numerical SCFT 



VIII. DISCUSSION 
A. Analytical mean-field approach 

The analytical mean-field theory is restricted in its va- 
lidity, whenever a lamellar phase is addressed, to small 
lamellar order-parameter amplitudes or to the proximity 
of a continuous microphase transition. In any case, it is 



The SCFT method invokes the mean-field approxi- 
mation, too, but avoids the assumption of small order- 
parameter amplitudes and the single-harmonic approx- 
imation for the lamellar phase. Thus, it provides ap- 
propriate mean-field predictions for large regions of the 
phase diagram, but due to numerical problems fails as the 
multicritical point is approached and both wave numbers 
and free-energy differences decrease. Moreover, numer- 
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ical SCFT is restricted to a small number of different 
components, and consequently allows us to address ran- 
dom copolymers with a small number of blocks Q only, 
which led to the choice Q — 3 in this study. The SCFT 
calculation for random continuous block copolymers with 
Q = 3 locates the entire three-phase region in the half- 
plane A > A c of the X~x diagram. 



C. Combining the results 

Beyond the mean-field approximation, the analytical 
approach and SCFT have different additional limitations, 
such that their results for the location of three-phase 
boundaries are complementary: The analytical approach 
assumes the lamellar order-parameter amplitudes to be 
small, which is accurate in the vicinity of the multicritical 
point. In this region, however, also the free-energy differ- 
ences between competing states (global lamellae, three- 
phase coexistence, two homogeneous phases) become mi- 
nuscule (cf. the inset of fig. [7]), which poses numerical 
difficulties for the SCFT calculations. Hence, there is 
no regime where both approaches are simultaneously re- 
liable, and a direct comparison is difficult. In the inset 
of fig. [14J we try to combine their results for the phase 
diagram of random continuous triblocks to one picture. 
The predictions for the cloud points of the homogeneous 
phases obtained by SCFT (dotted) and by the analytical 
method (dashed) match quite well, whereas the agree- 
ment for the cloud points of the lamellae is less satisfac- 
tory. Numerical SCFT results for these points (solid tri- 
angles) do not extend below A = —0.025 due to the men- 
tioned subtle free-energy differences in this region which 
control the phase behavior. The thin dotted line has not 
been computed, but marks our tentative extrapolation 
of SCFT data toward the multicritical point, based on 
the slope of the lamellar cloud line determined with the 
analytical theory (dot-dashed) in the part that is in qual- 
itative accordan ce ( thick). The analytical prediction for 
this line (cf. fig. |TTj ) is enhanced relative to the rougher 
description presented in fig. [7} cf. sec. VC Still, owing to 
the delicate free-energy balance, the shape of this bound- 
ary line is bound to be more sensitive to the approxima- 
tion of small lamellar amplitudes in the theory than that 
of the other three-phase boundary, at which the lamellar 
phase is the shadow and all amplitudes are smaller. 



we discover the coexisting phases to select sequences that 
match their morphology. Upon entering the three-phase 
region, the incipient shadow phase emerges with van- 
ishingly small volume fraction and with a sequence dis- 
tribution that already differs from the A-distribution of 
the cloud phase. Fractionation demixes the initial ran- 
dom (here Markovian) distribution into sequence classes 
(following the analytical approach, progressively), a sep- 
aration mechanism which might prove useful to isolate 
wanted species in polymer blends. 

Our analysis has been restricted to mean-field theory. 
For the macroscopic phase separation of the disordered 
state at A > A c , the critical region (x~ Xh)j within which 
the mean-field approximation fails, has been estimated 
with the help of a Ginzburg criterion [17]. The latter 
yields a Ginzburg number Gi oc Q 2 /M, i.e., the critical 
region does not shrink simply with chain length QM, in 
contrast to naive expectation. For fixed Q, such as con- 
sidered here, the mean-field predictions are correct in the 
limit of large M. The transition from the disordered to a 
global microphase-separated state (A < A c ), is expected 
to be weakly first-order due to fluctuations [HI [T3] . For 
the transition lines to three-phase coexistence and the 
multicritical point at A = A c , the effects of fluctuations 
remain to be explored. Whereas for simpler phase dia- 
grams, it has been shown that the Lifshitz point at A = A c 
is destroyed by fluctuations, the situation here is more 
complicated due to the fact that four phase states meet 
in the multicritical point. 

Phase coexistence enabled by component selection 
might be of interest for various other multi-component 
systems, cf., e.g., refs. [55] E3] - Specifically for polydis- 
perse copolymers, sequence fractionation can be general- 
ized starting from the case considered here: A straight- 
forward extension is to study random block copolymers 
asymmetric in global A-/B-content, which apart from 
the lamellar state, display other structured ordered mor- 
phologies, such as spheres on a bec lattice or hexagonally 
arranged cylinders [23J [35] . Other generalizations would 
include copolymers either with an arbitrary number of 
blocks or built from more than two segment types. Frac- 
tionation may also give rise to structured phases beyond 
the ordered microphases. Particularly promising in this 
context are random copolymers with many blocks, which 
might display frozen, random structures in coexistence 
with macroscopically phase-separated states. 



IX. CONCLUSIONS AND OUTLOOK 



The analytical method and the numerical SCFT con- 
stitute complementary approaches, which both have their 
virtues and together provide a comprehensive mean-field 
picture of the complex phase behavior of random triblock 
copolymers. With both methods, we consistently reveal 
an extended three-phase coexistence region of macro- 
scopic and microscopic phase separation in random tri- 
block copolymers, as suggested by simulations [TS]. Also, 
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APPENDIX 



Appendix A: Gaussian-chain averages 



Equation (22 1 and vertices of the expansion eq. (23) 



contain rt-point correlations of the Gaussian-chain mea- 
sure: 

^exp|-i^fc r r(s r )|\ = (Al) 

5j2 t fc,-,0 cxp < E l S r - Sr 'l fcr ' fe r/ \ 
[i<i' ) 

(derived for continuous chains in [36] , appendix B). 
Appendix B: Vertex functions and moments 



In eqs. (22 1 and (23 1, we also introduced the following 



functions: The discrete Debye function 

L 



D(L,k 2 ):= ]T (e 



si,s 2 = l 
L 



) -i(fci-7-(si)+fc2-r , (s 2 )) 



(Bl) 



_L(l + e~ k ) 2c- fc (l-c- Lk ) 
l-e- k2 ~ (1 -c- fe2 ) 2 ' 

and the structure factors, for individual sequences: 

L 



si,s 2 = l 
L 

= E 1^Ms2)e~ k2 ^- s ^, (B2) 



si,s 2 = l 



SW(fe l5 fe 2 ):= E ?v(*i)? v (*2)(e- i ^*'- r W 

Sl,S2,S3 = l 

(B3) 

S^(kx,k 2 ,k 3 ):= (B4) 



(B5) 



1 \s2 — Sl\—k 2 \si—S3\ 



Again, the length scale is the effective segment length b, 
k 2 := b 2 k 2 /(2d), with k the physical wave number. 



For the global A-distribution, the type correlation of 
two monomers on the same chain, whose block numbers 
differ by A(3(si, s 2 ) G {0, . . . , Q — 1}, can be calculated 
directly via the transition matrix M 

[q(si)q(s2)] x '■= Ypu{^)qu{si)qv{s 2 ) 

= A |A0( S1 , S2 )|_ (gg) 

Summing over all monomer pairs gives the second-order 
moment (cf. eq. (21)) for a A-distribution: 



™2(A) := ^ E fe( s iM s 2)] A 

si,s 2 = l 

0=3 1 2A (2-3A + A 3 ) 
3 + 9 (1-A) 2 ' 



(B7) 



Inserting eq. ( B6 1 into eq. ( 24 ) and performing the sum 



over all pairs yields the expression eq. ( 27 1 for the second- 
order structure factor S (fc 2 ) in a A-distribution. We 
abstain from presenting within this paper our compu- 
tations of the structure factors, eqs. (B2)-(B5), of a A- 



distribution for general Q (the expression eq. (27) had 
been given earlier in |37j). and of individual sequences 
for Q = 3. Obtaining the lengthy expressions for the 
fourth-order structure factors requires extended sorting 
of the multiple sums' terms due to combinatorics. 



Appendix C: Macroscopic phase separation 

Within the 'crushed polymer picture', we derive for the 
free energy of coexisting homogeneous phases a closed 
expression that is not limited to small order-parameter 
amplitudes. Here, each chain reduces to one structure- 
less particle with an A excess cjj equal to the average 
over all segments on that chain. Again with a field-based 
approach, the calculation of the free-energy functional is 



analogous to that in sec. Ill A but simpler, since confor- 
mational averages are obsolete for only one position Tj 
per chain. (For a replica-based derivation see [37] i the 
results prove to agree with Flory-Huggins theory [BJ.) 

For Q-block copolymers, it is sufficient to distinguish 
(Q + 1) components according to their A excess: 



qi 



21 -Q 



-qq-h I G {0,1,..., Q} 



(CI) 



In the case of symmetric triblock copolymers, the four 
component probabilities pi are related to the sequence 
probabilities defined in eq. ( 29 1 via 



Po=P3 = 



Pi 



Pi=Pi = 



P2 +P3 



1 



Pi 



2 ' r r 2 2 

With the coarse-grained component densities 



(C2) 



N 



ei {r) = L^8 5 . A 8(r-r j ), I e {0, 1, . . . , Q} , (C3) 

3 = 1 
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the total and A excess densities are 



and the partition function to calculate is 



N 



i=i 
exp 



oTVj 



(C5) 



^/^(x(-(r)) a -2«( ff (r)) a ) 



Introduction of additional fields, similarly as in eq. (13 1, 



and elimination of the original fields at the saddle point, 
yields the effective Hamiltonian per chain 



(C6) 



^/d^( x (^)) 2 -J(^)) 2 )-E^^ 
with the single- component partition functions 

*l := ~ J d d x exp j A ^f(x) - | . (C7) 

The general ansatz of JsT < (Q + 1) homogeneous phases, 



w(x) = = £ 



•,(fc) 



r(aj) = 7-W = ^}J« W ; , 



with volume fractions u^ 5 ) := |K| |/V gives 

_ _ (fe) f£( X afW-2o>W)\ 

z l=Z^ v W expj — — J \. (C8) 

We optimize ft with respect to the {v^ k \ uj( k \ r'*)}, 
with Lagrange multipliers Ai, As, A3 for the con- 
straints of number conservation, y~], — 1, constant 
global density, kqq — ^2 k v^Lb^ k \ and A excess, here 
v^t^ = 0. Solving the equilibrium conditions for 
nearly incompressible density conjugates, 



KQq 



1 + C^k" 1 + O («r 2 ) , 



(C9) 



we find expressions for the density-conjugate differences 

&(*>)- x (f( fc )) 2 - (f( fe ')) 2 



lim 

K-S-OO 



4 ^ 



(CIO) 



quadratic in the A-excess conjugates, such as in eq. (22 1 



Finally, we arrive at a self-consistent set of equations for 



the volume fractions and the field values, consisting of 
the constraints and 



go 



LX 
4 



f(fe) 
Qo 



go 



with the component partition functions (Cll) 



u fe ' exp 




and the {£>j } determined implicitly. 

For symmetric random triblock copolymers, the ansatz 
from eq. (41 1 of two homogeneous phases yields the am- 



plitude fh given in eq. (44) and, choosing Ai = —1, the 



free-energy density /h from eq. ( 43 1 



Appendix D: Equilibrium conditions at A > A c 

The parameter vector \y^ 2 \ n 2 , n 3 ] of the fractionated 
lamellar phase at A > A c must be determined as a zero of 
the gradient vector V^ (2) ni> n3 \/frac with components 



(1) 



f (2) _ Al) , "2 - P2 + n 3 - ps df± 

, , 1„ (l -P2 ~P3 - w (2) (1 - n 2 - n 3 )) 

f ) ?y In 



In 



(p y - v^n u ) (1 - n 2 - n 3 ) 
(l-n 2 -n 3 )(l-i;( 2 )) 
1 - P2 ~ Ps ~ v(2) (1 - n 2 - n 3 ) ' 



n 2 (l - P2 - P3 - « (2) (1 - n 2 - n 3 )) 



(Dl) 



In 



(p 2 - v( 2 )n 2 ) (1 — 712 — ^3) 



| 1-^)9/^ 
i9n 3 u( 2 ) 9ri2 

n 3 (I-P2- P3 - v {2) (1 - n 2 - n 3 )) 



In 



(p3 - v( 2 hi 3 ) (1 - n 2 - n 3 ) 



Here, p„ denote the constant A-defined concentrations 
p v (X), whereas n„ are the variable concentrations in the 
fractionated state. The following relations between the 
partial derivatives of f^\v^ 2 \ n 2l n 3 ) were inserted: 



dn 3 dn 2 
Ai) 



and 



dfb = n 2 -p 2 + n 3 - p 3 df^ 

dv (2) ~ v (2) _ v (2)} dn2 



(1) 



(D2a) 
(D2b) 
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The fact that depends on t he c oncentration (n 2 +n 3 ) 
only, simplifies the eq. system (Dl I for v^ 2 ' = const (e.g., 
in computing the three-phase boundaries; cf. sec. IV D I. 



From eqs. (|24j), (|39j), (|43l, and (|50j) one reads off the 
derivatives of ffp and fffl as functions of u^ 2 ), ri2, 11.3: 



9A 2) _ j 2 X (S„ (fc 2 pt ) - £> (3M, fe2 pt )) 



9n„ 



X S(fc 2 pt )-2L 



(D3) 



16 9m 



(2) 



81 «T 



,(2) 



/i 2) , ^ = 2,3, 



where A; 2 pt , an d TO i 2 '' are functions of n 2 , ^3, and 



(1) 



v (2) 



9n 2 



1 - «(a) 



In < 4 cosh 



6^0 



(i)' 



3>, (D4) 



with the amplitude ft /go (eq. (44 1 with pi — Uj 1 ) ex- 
pressed as a function of v^ 2 \ n 2 , 713 via eq. (50). 



At a given point (A, L\) of the phase space, the allowed 
domain V for the variables v^ 2 \ n 2 , is 



V : = 



,(2) 



n 3 e 



e [0,1], n 2 £ 



0, min ( 1, 



P2 



0, min ( 1 — ri2 



P3 

u( 2 ) 



(D5) 



fcopt(^2,n 3 ) > and Lx m (n 2 ,n 3 ) < L\ 



Appendix E: Numerical solution of the equilibrium 
conditions for the fractionation free energy 



In order to locate the zeros of the system (Dl ), which 



correspond to a minimum of the fractionation free energy, 
we employ a Newton-type algorithm using the following 
steps (exemplified for the fractionated lamellar phase): 



1. At a given (A,x), guess start parameter vector 

(2) 1 ^ 

xq := Vq ,n 2 , 0,^3,0 (A,x)- The sensitivity re- 
garding the start vector impedes completely autom- 
atized scans in the X-\ plane. 



2. Iteratively, apply Newton scheme 

X! = Xq - H^ 1 (x ) .V/frac (x ) 



(El) 



with H the Hessian of the system (Dl ) 



3. Stop if either the desired relative precision e := 
^I^J^ or a given maximal number of iterations 
has been reached. In the latter case, and if H gets 
singular during the iteration, restart from step [T] 

4. To ensure that /f rac (sci) is a minimum, check H for 
positive definitencss, i.e. calculate its eigenvalues. 

5. From the minimum concentrations n 2 , n 3 , calcu- 
late fc op t(ri2) ^3) and Lx m (n 2 , 71.3) (the continuous 
microphase transition that would occur in an inde- 
pendent, global sequence distribution equal to the 
fractionated one) and ensure the result vector to 



comply with eq. (D5I 



Convergence, especially while approaching the multi- 
critical point (\ c ,Lxc), can be achieved only for start 
vectors very close to the actual solution. Therefore, pro- 



ceeding on a three-phase boundary line (see section IV D ) 
toward (X c ,Lx c ), we use the solution at one value of A 
as the start vector for the adjacent A. The resolution for 
A is chosen between 5 x 10~ 4 far from A c and 10 -5 near 
A c , and between 10 -3 and 10~ 4 for L\. In the vicinity 
of (A c , Lxc), entries of the start vector have to be even 
closer to the actual solution and are obtained by extrapo- 
lating solutions on the boundary line. Finally, the result 
vector is calculated with a relative precision e = 10~ 12 
of its modulus. Uniqueness of solutions of the nonlinear 



equation systems mentioned in sees. IV B IV D cannot 
be proven rigorously. However, we are sure not to miss 
transition lines to three-phase coexistence at lower L\, 
since at each A, we start to scan the domain of definition 



eq. (D5) with the A-defined concentrations. 
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